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1.  Introduction  and  definitions 

One  of  the  main  goals  of  every  science  is  the  endeavoring  to  arrange  and 
collate  the  numerous  individual  observations  and  details  which  present  themselves, 
in  order  that  they  become  part  of  one  comprehensive  picture.**  Such  an  objective 
may  be  achieved  through  the  use  of  data  reduction  which  is  an  essential  part 
of  chemical  and  physical  ciata  analysis.  In  the  past  the  task  of  data  reduction 
involved  mainly  curve  fitting  procedures  which  have  recently  become  more 
efficient  due  to  increased  availability  of  computers.  Topological  approaches 
which  consider  additivities  of  molecular  properties  have  been  almost  ignored 
in  the  past.  However  some  recognition  of  the  role  of  Rraph  theory  became 
apparent  soon  after  the  (important)  work  of  Smolensk ii^  and  (to  a  lesser  degree) 
the  work  of  Gordon  and  Kennedy^  who  formalized  their  schemes  of  general 
approaches  to  molecular  additivities.  Revival  of  interest  in  chemical  graph  theory 
is  probably  due  to  Balaban  via  his  (editing)  of  a  book4  which  provided  a  source 
of  problems  and  numerous  other  papers.  Several  counting  polynomials  were  defined 
which  act  as  descriptors  of  the  connectivity  pattern  of  various  types  of  molecular 
graphs.  The  latters  are  well-defined  mathematical  objects  which  describe  the 
bonding  relation  between  atoms  in  the  original  molecules.  Such  counting 
polynomials  are  then  a  combinatorial  form  of  data  reduction.  Equation  1  describes 
the  general  form  of  a  graph- theoretical  counting  polynomial^,  f(C;x),  for  an 
arbitrary  graph  G,  containing  n  vertices: 

M 

F(G;x)  =  £  p  0(G;k)  xf(k'n)  (1) 

k=0 

where  P  is  either  1  or  (-1)k,  0(G;k)  enumerates^certain  selected  graph-invariants 
in  G  taken  in  k  independent  tuples  (i.e.  no  two  of  them  are  adjacent)  and  M  is 
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the  maximal  value  of  k.  For  the  purpose  of  this  paper  we  consider  three 
polynomials  for  all  of  which  p  •  1  and  f(k,n)  ■  n.  These  are: 

i)  The  counting  polynomial  (of  Hosoya)®,  H(G;x)  for  which  6(G;k)  «  p(G;k)  * 
the  number  of  k-matchings  in  G. 

ii)  The  sextet  polynomial  (of  Hosoya  and  Yamaguchi)?,  o(B;x)  of  which  6(G;k) 
-  r(B;k)  *  the  number  of  selections  of  k  mutually  resonant  but  nonadjacent 
aromatic  sextet  in  a  benzenoid  system  ( =  polyhex  graph®),  B. 

iii)  The  independence  polynomial  (of  Gutman)^,  u(G;x),  for  which  6(G;k)  *  0(G;k) 
=  the  number  of  selections  of  k  independent  vertices  of  C. 

2.  The  Caterpillar  Tree^  the  Polyhex  Graph”'*'  and  the  Clar  Craph9»”l2 


In  a  given  Kekute  structure  an  (aromatic)  sextet  is  defined  as  a  set  of  three 
double  bonds  circularly  conjugated  so  that  no  two  sextets  can  have  a  common 
bond.  When  two  rings  (hexagons)  can  have  aromatic  sextets  so  that  all  other 
carbon  atoms  are  spanned  by  a  sextet  or  by  double  bond,  such  two  rings  are  called 
resonant. 9,12  The  individual  hexagons  of  a  given  benzenoid  system  may  or  may 
not  be  resonant^.  Information  on  such  resonance  relations  (among  the  individual 
hexagons)  are  best  "reduced”  (i.e.  described)  using  either  a  caterpillar  tree^ 
or  a  Clar  graph. ^  To  illustrate  this  fundamental  relation”'^  we  consider  the  three 
graphs  shown  below: 
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The  hexagons  of  B  may  be  grouped  into  three  subsets,  viz.,  (1,2,3);  { 4, 5,6, 7,8) ; 
{9,10,11,12).  One  observes  that  no  two  hexagons  are  resonant  if  they  belon 


to  one  subset.  Similarly  the  edges  of  T  may  be  subdivided  into  analgous  three 
subsets  observing  that  no  two  edges  of  the  same  subset  can  be  adjacent.  A  similar 
description  obtains  for  the  vertices  of  graphy^.  In  fact  bijective  (i.e.  one-to-one, 
onto)14  mapping  might  be  defined  between  the  hexagons  of  B,  the  edges  of  T 
and  the  vertices  of  A:  the  latter  is  called  the  Clar  graph9,12  of  b.  It  is  easy 
to  verify  the  following  identity; 


o  (B;x)  -  H(T;x)  «  w<A;x)  *  1  +  12X  ♦  45X2  ♦  53X3 


A  caterpillar  tree  such  as  T  whose  counting  polynomial  is  identical  to  the  sextet 


lynomial  of  a  benzenoid  hydrocarbon  B  will  be  called  an  equivalent  caterpillar 


to  the  benzenoid  hydrocarbon,  the  set  of  graphs  {  B,T,A)  shown  above  will  be 
called  a  set  of  equivalent  Rraphs.  Gutman^3  demonstrated  that  for  every 
nonbranched"*^  benzenoid  hydrocarbon  there  is  an  "equivalent"  caterpillar  tree. 

3.  Conjecture 

The  physical  properties  of  a  benzenoid  hydrocarbon  which  depend  on  its  topological 
structure  may  be  predicted  from  graph-theor<ftical  properties  of  its  equivalent 


4.  Results  and  tests  of  Conjectures 

Fig.  1  shows  benzenoid  hydrocarbons  and  their  equivalent  caterpillar  trees  studied  in  this 
4.a  Graph  theoretical  (combinatorial)  properties 


Three  graph -theoretical  properties  of  benzenoid  hydrocarbons,  viz., 
self*avoiding  paths^,  conjugated  circuits^  and  the  number  of  Kekute  structures 
(i.e.  number  of  perfect  matchings)  are  correlated  with  the  molecular  connectivity 
indices, 21  x 's  of  the  equivalent  caterpillar  (or  pseudocaterpillar)  trees.  Fig. 
2  shows  the  types  of  linear  correlations  obtained.  In  all  cases  correlation 
coefficients  >  0.999  were  obtained. 


4b.  Physical  properties 

Two  physical  properties  viz.,  electronic  absorption  spectra 22,  and  heats 
of  atomizations^  of  several  homologous  series  of  benzenoid  hydrocarbons  are 
linearly  correlated  with  x  CD's;  the  molecular  connectivity  indices^  of  the 
corresponding  equivalent  trees.  Details  of  data  are  available  oh  request ^  but 
as  an  illustraton  Table  1  is  included  which  contains  6  and  p  bands  of  the  electronic 
spectra  of  two  families  of  hydrocarbons  of  Fig.  3.  A  separate  manuscript 

is  submitted  elsewhere^4  dealing  with  electronic  absorption  spectra  of  27  arene 
systems  all  linearly  correlated  withJ^(T)  values.  The  same  type  of  data  reduction 
is  obtainable^  when  molecular  properties  (topological  or  physical)  of  benzenoid 
hydrocarbons  are  studied  in  terms  of  connectivity  functions  of  their  equivalent 
Clar  graphs.25  Fig.  4  is  a  plot  of  heats  of  atomizations  of  some  benzenoid 
hydrocarbons  versus  connectivity  indices  of  the  corresponding  caterpillar  trees. 


(. 

0 

6.  Discussion 


Caterpillar  trees  and  Clar  graphs  seem  to  be  quite  promising  "storage 
devices"  which  preserve  many  properties  of  (the  much  larger)  equivalent  polyhex 
graphs.  In  contrast  to  the  Balaban-Harary  dualist  graphs^  they  do  not  have 
a  geometric  element  and  as  such,  connectivity  alone  best  describes  them.  It 
is  amazing  that  some  of  the  studied  properties  such  as  the  total  number  of  paths^9 
in  a  polyhex  graph  (Fig.  2)  take  nearly  30  minutes  on  a  Corona  computer, 27 
yet  can  be  retrieved  by  very  simple  hand  calculation  of  the  connectivity  index 
of  a  much  simpler  graph!  Table  2  lists  percent  retrievals  of  the  total  number 
of  self-avoiding  paths  for  a  series  of  nonbranched  benzenoid  hydrocarbons  shown 
in  Fig.  1.  All  the  above  mentioned  correlations  imply  classification  of  molecules 
into  families  (i.e.  homologous  series).  The  concept  of  classification  is  well-known 
in  science:  the  periodic  table  of  the  elements  and  recently  table  of  alkanes 
introduced  by  Randid  and  Wilkins^  and  more  recently  table  -  of  benzenoid 
hydrocarbons  of  Dias^^  are  all  non-mathematical  forms  of  set  theory.  However 
in  constrast  to  the  latter  two  tables,  which  result  from  graph  theoretical 
considerations,  the  "old"  periodic  table  is  a  result  of  group  thecfry  as  it  follows 
from  symmetry  properties  of  atoms. 
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Table  1 

Electron  absorption  spectra  of  the  6  bands,  Xg  ,  of  the  homologous  series  of 
benzenoid  hydrocarbons  shown  in  Fig.  3.  X  (T/s  are  molecular  connectivity  indices 
of  the  equivalent  pseudocaterpillar  trees. 


Hydrocarbon* 

X(T) 

(A) 

*  p( 

A 

3.914 

2570 

B 

4.684 

2870 

C 

5.293 

3130 

D 

5.812 

3430 

1 

1.000 

2068 

2 

1.414 

2850 

3 

1.732 

3745 

4 

2.000 

4710 

5 

2.236 

5755 

6 

•c.f.  Fig.  3 

2.449 

6930 

Table  2 
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Percent  retrieval  of  the  total  number  of  self-avoiding  paths  for  the  series  of 
nonbranched  benzenoid  hydrocarbons  shown  in  Fig.  2.  The  regression  line  is  Y 
*  4.374  ♦  0.484X  with  a  correlation  coefficient  =  1.0000. 
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Polyhex  graphs  of  benzenoid  hydrocarbons  studied  in  this  work  together  with 
the  molecular  graphs  of  their  corresponding  caterpillar  trees. 


Fi&^2 

Combinatorial  properties  of  benzenoid  hydrocarbons  studied:  K(B),  Ep(B)  and 
I  Rtot*®*  respectively  indicate  number  of  Kekul£  structures,  total  number  of 
self-avoiding  paths  and  total  number  of  conjugated  circuits  of  a  benzenoid  system, 
B,  plotted  versus  the  connectivity  index  of  the  corresponding  caterpillar  tree, 
X(T)  raised  to  various  powers.  (The  hydrocarbons  and  their  caterpillars  are  shown 
in  Fig.  1). 

PiR-  3 

Electronic  absorption  spectra  of  benzenoid  hydrocarbons:  Points  1-6  indicate 
benzene,  naphthalene,  anthracene,  tetracene,  pentacene  and  hexacene;  the  line 
passing  through  these  points  is  a  plot  of  their  absorption  para  band  (in  In  units) 
vs.  the  connectivity  index  of  the  corresponding  caterpillar  trees,  *(T).  The  line 
passing  through  points  A-D  is  a  plot  of  the  beta  band  vs.  X  (T)  where  the 
hydrocarbons  and  their  trees  are  shown  in  Fig.  1. 

PiR.  4 

A  plot  of  the  natural  logarithms  of  the  heats  of  atomizations  of  some  benzenoid 
hydrocarbons.  In  Ht  At(B),  against  the  connectivity  indices,  X(T/s,  of  the 
corresponding  caterpillar  trees  shown  in  Fig.  1. 
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